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ABSTRACT 

We have observed the Sunyaev-Zel'dovich effect in a sample of five moderate-redshift 
clusters with the Ryle Telescope, and used them in conjunction with X-ray imaging and spec- 
tral data from ROSAT and ASCA to measure the Hubble constant. This sample was chosen 
with a strict X-ray flux limit using both the BCS and NORAS cluster catalogues to be well 
above the surface-brightness limit of the ROSAT All-Sky Survey, and hence to be unbiased 
with respect to the orientation of the cluster This controls the major potential systematic ef- 
fect in the SZ/X-ray method of measureing Hq. Taking the weighted geometric mean of the 
results and including the main sources of random error, namely the noise in the SZ measure- 
ment, the uncertainty in the X-ray temperatures and the unknown ellipticity of the clusters, we 
find//o = 59+7kms"'Mpc"' assuming a standard CDM model with i^M = 10,i2A = 0.0, or 
Ho 



:65;^kms-iMpc-^ if iijvi = 0.3,£2a = 0.7. 
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1 INTRODUCTION 

It ha s long been recognised that th e Sunyaev-Zel'dovich (SZ) ef- 
fect (Sunyaev & Zel'dovich 1972) in clusters of galaxies, along 



with X-ray measurements, can p rovide a method of measuring dis- 



Cavaliere et al. 1975 



tances on cosmological scales (junn 1978; Silk &. White 1978 
Birkinshaw 1979| ). The method takes advan 



tage of the differing dependences on the cluster electron density 
We and temperature Tg of the X-ray bremsstrahlung n^.T^^'^) and 
SZ effect (oc «e7i) to measure the physical size of the cluster, and 
hence its angular-diameter distance. The dependence of the derived 
value of Hq on the important quantities is 

\ATszJ ^11 

where AJsz is the SZ brightness decrement, 5x is the X-ray surface 
brightness, 6 is the apparent angular size of the cluster, and is 
the ratio of the linear sizes of the cluster perpendicular and parallel 
to the line of sight. 

The systematic effects in this method of measuring Hq are 



quite different from those in other methods, either in the traditional 
distance ladder or in other direct physical methods such as gravita- 
tional tensing, and are mostly due to the difficulty of modelling the 
cluster gas that causes both the SZ effect and the X-ray emission. 
Recent work has shown that the small-scale features in the gas dis- 
tribution such as cooling flows or clumping contribute smaller bi- 
ases than might have been thought, since the effects on the observed 
quantities entering H q tend to cancel each other even for a single 
cluster (eg, Grainger ( 2001 ), Maggi et al.). However, the large-scale 



distribution of the gas poses a problem, since the line-of-sight depth 
through the cluster cannot be measured directly, and has to be as- 
sumed to be the same as the size of the cluster in the plane of the 
sky. This asumption will be violated if clusters are chosen prefer- 
entially because of high surface brightness, since this will favour 
clusters that are extended along the line of sight. This will bias low 
the derived value of Hq. We therefore need a cluster sample that is 
selected on total flux rather than surface brightness, well above the 
surface brightness limit of the survey. Although several individual 
Hj) estimates from the SZ/X-ray method have been p ublished (eg 
(jHughes & Birkinshaw 1998|; ^olzapfel et al. 1997; Reese et al. 
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2000)), properl y selected cluster samples are onl y recently begin- 
ning to emerge ( Mason, Myers & Readhead 200 1 ). In this paper we 



derive a sample of moderate-redshift X-ray-selected clusters that is 
carefully controlled to be free of any orientation bias, and use it to 
calculate //y using the SZ data from the Ryle Telescope (RT). 

Unless stated otherwise, we assume a standard cold dark mat- 
ter cosmology (SCDM, Sly^ = 1.0, Q.^ = 0.0). We also give our 
results assuming a ACDM cosmology (flM = 0.3, Slj^ = 0.7). 



2 X-RAY CLUSTER SELECTION 

The most complete X-ray cluster samples suitable for our purpose 
are selected from the ROSAT All-Sky Survey (RASS). We use 
two s_uclr^am£leSj_^he^ri^h£^^ et al. 

1998), and Northern Rosat All-Sky Survey (NORAS) (B ohringer et 
al. 2000). We began selection of our SZ sample from BCS while it 
was being compiled; we use NORAS as a consistency check. Both 
surveys are based on cluster candidates found in the first processing 
(SASS) of the RASS data, with further tests to check they are in- 
deed clusters. A problem, particularly for clusters at high redshift, 
is the risk of incompleteness for compact clusters that may not ap- 
pear extended enough to be recognised as such. BCS attempts to 
overcome this problem by adding known Abell and Zwicky clusters 
to the sample. However, comparing the BCS and NORAS samples 
after applying our RT selection criteria, it turns out that our sam- 
ple is not affected by differences between the BCS and NORAS 
secondary selections. 

The flux limit of NORAS is 3 x 10"'^ erg s"' cm^^ at 90% 
completeness; that of BCS is 4.4 x 10^'^ erg s^' cm^^ to the same 
completeness. But we limit our sample to even higher fluxes, of 
greater than 5 x 10^'^ erg s^' cm^^ (see Figure |^. Both this and 
a restriction to redshifts z < 0.3 help reduce the risk of missing 
compact clusters. In accordance with the requirements for the Ryle 
Telescope we restrict our sample to the northern sky at declina- 
tions 5 > 20°, to a redshift range of z > 0.14, and luminosities of 
L.v > 8 x 10'*'* erg s"' . Table|l] lists the clusters from BCS and from 
NORAS that meet these criteria, and also indicates those which 
have radio sources that would make SZ measurement difficult or 
impossible. Where the surveys disagree on whether a cluster should 
be included on grounds of flux, we cautiously select only those 
clusters that meet our criteria in both surveys: this excludes three 
clusters above our flux limits in BCS and not in NORAS, and one 
in NORAS but not in BCS. This is consistent with the errors in the 
flux estimates in each survey. Two BCS-selected clusters are not in 
NORAS, consistent with the inclusion of optically selected clusters 
in BCS that are not in NORAS due to lack of extent. Of the eleven 
clusters satisfying our X-ray criteria, six have radio sources in the 
cluster or within the field of view of the RT that are too bright at 
15 GHz to allow SZ observation with the RT. We have observed 
the remaining five, A697, A773, A1413, A1914 and A2218 with 
the RT and obtained detections for all of them; the results are dis- 
cussed below. Since the presence of radio sources cannot be cor- 
related with cluster orientation, this subsample is effectively ran- 
domly chosen, and therefore retains the orientation-independence 
of the whole sample. Indeed, we have checked that the our subsam- 
ple is not correlated to X-ray properties such as flux, extent like- 
lihood and measured extent. It therefore appears unlikely that any 
significant bias, larger than the random error, is present in our aver- 
aged Hq result. Clearly a larger sample would improve matters: the 
small size of our sample highlights the difficulty of selecting rig- 



orously at moderately high redshift, even using the best large-area 
X-ray surveys available. 



3 SZ OBSERVATIONS AND Ho CALCULATION 

Details of our observ ational and mod elling m ethods are given in 
Grainge et al. (2001a), Grainge et al. (2001b) and Grainger et al. 
(2001). We briefly summarise them here. 

In each case we fit the ROSAT PSPC or HRI image using an 
ellipsoidal P-model, with P, the core radii in the plane of the sky, 9i 
and 62, the position angle of the major axis, the position on the sky, 
and the central electron density no all as free parameters. The core 
radius in the line of sight is set to the geometric mean of the other 
two. A downhill simplex method is used to maximise the likelihood 
of the data given the model, using Poisson statistics. For the PSPC 
images we use the hard band (0.5-2 keV) in order to minimise the 
effect of Galactic absorption. 

We remove contaminating radio sources by a maximum like- 
lihood fit to the visibility data, simultaneously fitting the SZ decre- 
ment (using the X-ray-determined model, with a variable ampli- 
tude) and the sources. We include trial sources at all positions sug- 
gested either from the l ong baseline RT data or from other obser- 
vations such as NVSS (Z^ondon et al. 1998), and allow their po- 



sitions and fluxes to vary in the fit. This allows us to remove the 
effects of sources whose existence is known from lower-frequency 
observations but whose brightness would not warrant a significant 
detection from the RT data alone. 

Finally we take the source-subtracted data and compare them 
with mock SZ data derived from the X-ray model, taking into ac- 
count the full response of the RT. We vary the normalisation of 
the mock data (equivalent to varying Hq) and find the likelihood of 
the data, and then multiply by a prior that is uniform in log space 
(since Hq is a scale parameter) to find the peak and extent of the 
posterior probability distribution. We then apply corrections for the 
more exact relativistic_^om^^f^he_Com£toi^j^_^ara^^ 
nor & Lasenby 1998), and small time-dependant variations from 
the nominal RT flux calibration scale (determined from VLA mon- 
itoring of our primary flux calibrators, 3C48 and 3C286), as well 
as adding in (in quadrature) estimates of the other main sources of 
error. 

We now give details of the procedure for each of the five clus- 
ters. Two have b een published more fully elsewhere (A773, Saun- 
ders et al. ([200 1[) and A1413, Grainge et al. ([2001b|)) and are briefly 



summarised here; one ( A2218) is a re-analysis of previously pub- 
lished data (tones 19951); the other two (A697 and A 19 14) are new. 



3.1 A697 

A697 was observed with the ROSAT HRI on 1995 November 6 
with a live time of 28.1 ks. We also observed it using ASCA dur- 
ing 1996 April in order to determine its temperature. Analysing 
both the GIS and SIS data using standard XSPEC tools, we find 
a temperature of 8.0^q 4 keV (68% confidence), with a metallicity 
of 0.24 ±0.07 solar. Fitting to the J?OSAT HRI data, we find the 
parameters given in Table ^ The HRI image, model and residuals 
are shown in Figure ^ 

We observed A697 with the RT on 42 occasions between 
1994 August 6 and 1996 April 27. The map of all the data is 
shown in Figure ^ The bright source to the NW is extended, but 
can be successfully modelled as three point sources. Fitting si- 
multaneously to the sources and SZ decrement gives a good fit 
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with the source parameters listed in Table ^ Fitting the source- 
subtracted data to the model derived from the X-ray data (Fig- 
ure |3|), and applying the relativistic and flux-scale corrections, we 



and H() = 1 19^33 kms ' Mpc ' including the other sources of er- 



find//o = 52^^kms 'Mpc ' (errors from SZ fitting only). Includ' 



ing the other errors, dominated by the X-ray temperature (^i()%) 
and ellipticity (±14%), we find Hq = 52+] j kms"' Mpc"' . For a 
ACDM cosmology with Q.f^ = 0.3 and Q.\ = 0.7, this becomes 



//o = 59+}|kms 



"'Mpc"' 



3.2 A773 

Details of th e dete rmination of from A773 are given in Saun- 
ders et al. ( |200l| ). We used a ROSAT HRI image and a tem- 
perature measurement from ASCA in conjunction with our RT 

_[^kms"'Mpc"' forSCDM, 
'Mpc"' for ACDM. 



data to estimate a value of Hq = 77^ kms 



:85+^^kms- 



3.3 A1413 



Details for A1413 are given by Grainge et al. (2001bl, who also 
discuss the contributions to the Hq error budget in general. Us- 
ing a ROSAT PSPC image and ASCA temperature, we find Hq = 
for SCDM, Hq = 61+f^kms"' Mpc"' for 



57+^^ km s"' Mpc" 
ACDM 



3.4 A1914 

A1914 was observed with the ROSAT PSPC for 8.6 ks in 1992 
July. There are several pieces of evidence to suggest that A1914 
is undergoing a merger. The ROSAT image (see Figure ^ shows 
significant substructure in the core of the cluster, but no evidence 
of a cooling flow. The temperature, measured from an ASCA ob- 
servation, is rather high at 8.56 ±0.39 keV. The POSS optical im- 
age shows two distinct groups of galaxies with no single dominant 
galaxy, and finally there is a diffuse steep-spectrum radio source to 
the NE of the cluster centre that is detected in the WENSS (300 
MHz) ( |Rengelink et al. 1996| ) and NVSS (1.4 GHz) surveys. This 
raises the question of whether A1914 should be included in our 
sample, since it is possible that the cluster gas might be far from 
hydrostatic equilibrium and difficult to model. We cho ose to in- 
clude it in our sample since simulations (Grainger 2001) show that 
even including clusters with large internal kinetic energies in their 
gas does not bias the sample average Hq However, such clusters can 
have a large scatter about the mean value - this does appear to be 
the case with A1914 and also with A2218. 

We fitted a (3 model to the ROSAT PSPC image with the pa- 
rameters given in Table ^. The residuals (see Figure ^ were no- 
ticeably poorer than those from the fits to other clusters such as 
A1413; a Monte Carlo analysis shows that the hypothesis that the 
true structure is a p model and that the features seen result from 
Poisson noise can be rejected at the 4-a level. A1914 was observed 
with the RT on 19 occasions between 1996 March 6 and April 
16; the full resolution radio image is shown in Figure ^. There 
are clearly several radio sources near the cluster centre, some of 
which are also detected in a VLA image (Figure We used the 
VLA source positions and positions from a long-baseline RT im- 
age to estimate the fluxes of nine sources; the fluxes found and sub- 
tracted are shown in Table ^. Fitting the source-subtracted short- 
baseline data to the X-ray-based model (Figure ^) results in a value 
of //q = 1 19+3jkms"' Mpc"' with errors just from the SZ fitting. 



-38' 

ror. For ACDM this becomes Hi 



' 



129+^^ kms"' 



Mpc 



3.5 A2218 



A2218 was the first clus ter to be observed in the RT SZ pro- 
gramme ( Jones et al. 1993 ), and has also been extensively observed 



in many wa veban ds, with SZ detections in the radio (eg Birkinshaw 
& Hughes (1994) ) as w ell as strong and w eak len sing observations 
(eg Kneib et al. ( 199j) an d Squires et al. (1996)) and X-ray stud- 
ies (eg Markevitch (1997)). We h ave previou sly reported a mea- 
surement of Hq from this cluster (Jones 1995); we re-analyse our 



data here because we have improved our source-subtraction and 
modelling techniques, and because a new temperature measure- 
ment from ASCA is available that supersedes the previous GINGA 
measurement. 

Figure ^ shows the image made by combining 30 days' RT 
data taken in two array configurations between 1992 December 4 
and 1993 March 23. This map has simply been CLEANed with no 
source subtraction; nevertheless, even at this high resolution the SZ 
decrement is clearly visible, and an approximate calculation of the 
observed temperature decrement, using the CLEAN restoring beam 
size of 40" X 36", gives a maximum decrement of — 760± 150(UK. 
This agrees well with the fitted P-model central temperature of 
-797 |yK. 

We used the ROSATPSPC observation made on 1991 May 25, 
which has a live time of 42.47 ks. The X-ray model fitting parame- 
ters are shown in Table ^ As with A1914, the X-ray image contains 
substructure that is significantly different from a smooth ellipsoidal 
P model (see Figure p|). Preliminary source subtraction in the map 
plane had shown that, of the three bright sources evident in Fig- 
ure ^, the central one was point-like but the eastern and western 
sources had some additional extended flux. We therefore fitted a to- 
tal of six trial sources, modelling the western source as three points 
and the eastern source as two. This gave a good fit to the data with 
the source parameters listed in Table ^ Subtracting these sources 
and fitting to the X-ray m odel (using an ASCA-derived te mper- 
ature of 7.20 ±0.36 keV ( jVIushotzky & Loewenstein 1997| )), we 
find a best fit value of Hq = 3 1 +g km s" ' Mpc" ' , where these errors 
only include the SZ fitting. Including the other dominant sources 
of error (±10% from the temperature measurement and ±14% 



from the likely ellipticity), we find Hq = 31+g'^kms ' Mpc or 
Hq = 35l^' kms"' Mpc"' for ACDM. 



3.6 Combined result 

The values of Hq derived from each cluster are summarised in in Ta- 
ble ^ Since the main sources of error in each value are in the form 
of a multiplicative factor, we calculate the sample mean by taking a 
weighted average of log/Zy. The result is Hq = 591^ kms"' Mpc"' 
for SCDM or Hq = 65+^ kms"' Mpc"' for ACDM. Looking at the 
scatter in our individual values, we see that A697, A773 and A1413 
all agree with the mean within l-o. A1914 is 1.6 a above the mean 
and A2218 is 2.8 O below it. These are the two clusters for which 
there is independent evidence, from the lack of a good fit of the 
X-ray model and from the optical galaxy distribution, of signifi- 
cant dynamical activity. Since our error estimate does not include 
a component to take account of this effect (for lack of a means of 
quantifiying it), it is not suprising that the scatter appears anoma- 
lously large. Nevertheless we feel justified in keeping these two 
clusters in the sample based on the evidence from simulations that 
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dynamically active clusters do not bias the value of Hq - indeed, 
removing them does not significantly change the best estimate of 



4 CONCLUSIONS 

We have used ROSAT and ASCA X-ray data and Ryle Telescope 
SZ observations to measure the Hubble constant Hq from a sam- 
ple of five clusters, selected using the BCS and NORAS cluster 
catalogues to be free of orientation bias. The weighted geometric 
mean value of Hq from this sample is H^) — 59^7kms^' Mpc~' 
assuming a standard CDM model with = 1.0, £1a = 0.0, or 
Hq = 65l^ kms-i Mpc-i if Hm = 0.3, Ha = 0.7. 
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Cluster 


Flux density (jj}y) 


RA offset (arcsec) 


dec offset (arcsec) 


A697 


462 


-468 


193 




154 


-442 


156 




215 


-483 


226 




624 


28 


-236 




1 <1 


— ZOZ 


— lyj 




no 


1 /^Q 

— loy 


QQ 




155 


-141 


31 




173 


-108 


201 




103 


-118 


-52 




278 


-96 


-82 




234 


49 


66 




123 


-56 


210 




186 


13 


158 


A2218 


1438 


-189 


36 




1227 


-40 


117 




1541 


126 


94 




522 


132 


105 




543 


-196 


23 




240 


-207 


41 



Table 3. Sources subtracted from the radio data. The flux densities given 
are apparent, ie not corrected for the RT primary beam. Offsets are 
from the radio pointing centres: A697 08'>42'"57\7 + 36°.721'45"; A1914 
14''26"02\0 + 37°49'29"; A2218 14''26""02\0 + 37°49'29"(J2000). The 
errors on the source fluxes are ±30, ±45 and ±35 /Jy respectively. 



Cluster Ho (SCDM) Hq (ACDM) ATb (/jK) 



A697 


52+13 


59l|i 


1047 


A773 


77+19 
"-17 




737 


A1413 


57+23 
^'-16 


61+25 


863 


A1914 




129+51 


864 


A2218 


31^r 


35tl' 


797 



Sample 59^, 65+^ 



Table 4. Hq values for the five clusters in the sample, for both SCDM and 
ACDM cosmologies. Also given are the central temperature decrements us- 
ing the best normalisation of the fitted P model to the SZ data. 
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Cluster 


R.A. (J2000) 


Dec 


z 






Jx 


, MOD A 


B? 


N? 


R? 


A586 


113.091 


31.629 


0.1710 


9.1 


11.12 


7.94 


9.84 


* 


* 




A665 


127.739 


65.854 


0.1818 


11.8 


16.33 


11.23 


15.69 


* 


* 




A697 


130.741 


36 365 


0.2820 


5.0 


16.30 


5.77 


19.15 


* 


* 


* 


A773 


139.475 


51.716 


0.2170 


6.7 


13.08 


5.71 


12.11 


* 


* 


* 


Z2701 


148.198 


51.891 


0.2140 


5.6 


10.68 












A963 


154.255 


39.029 


0.2060 


5.9 


10.41 












A1413 


1 78 826 


23.408 


0.1427 


15.5 


13.28 


12.69 


10.91 


* 




* 


A1423 


179.342 


33.632 


0.2130 


5.3 


10.03 


3.73 


7.23 


* 




* 


A 1682 


196.739 


46.545 


0.2260 


5.3 


11.26 


4.07 


8.79 


* 






A1758 


203.179 


50.550 


0.2799 


3.6 


11.68 


5.60 


18.29 




* 


* 


A1763 


903 81 8 


40.996 


0.2279 


6.9 


14.93 


6.35 


13.85 


* 


* 




A1914 


216.509 


37.835 


0.1712 


15.0 


18.39 


12.90 


15.91 


* 


* 


* 


A2111 


234.924 


34.417 


0.2290 


5.0 


10.94 


3.91 


8.68 


* 






A2218 


248.970 


66.214 


0.1710 


7.5 


9.30 


7.16 


8.16 


* 


* 


* 


A2219 


250.094 


46.706 


0.2281 


9.5 


20.40 


11.18 


24.33 








RXJ1720 


260.037 


26.635 


0.1640 


14.3 


16.12 


10.84 


12.34 




* 




A2261 


260.615 


32.127 


0.2240 


8.7 


18.18 


9.83 


20.61 


* 


* 





Table 1. The selection of the cluster sample from the BCS and NORAS surveys. Clusters are required to have Lx > 8 x 10^'*ergs^' and /x > 5 x 
10~'^ergs~' cm~^ in both samples, and to be free of contaminating radio sources. The final three columns indicate which clusters are selected from BCS, 
NORAS, and radio data respectively. Eleven clusters meet the X-ray selection criteria, and five of these (shown in bold) are sufficiently free of radio sources 
to have their SZ effects measured. Units of luminosity Lx are 10'**ergs~' and of flux fx are 10~'^ergs~' cm~^. 



Cluster 


z 


P 


ei (") 


62 (") 


PA° 




Temperature (keV) 


A697 

A773 
A1413 
A1914 
A2218 


0.282 
0.217 
0.143 
0.1712 
0.171 


0.67 

0.64 
0.58 
0.69 
0.65 


46 

61 
56 
51 
69 


34 

44 
38 
46 
53 


101 

16 
1 

10 

-8 


7.9 X 10-3 

6.8 X 10-3 
9.6 X 10-3 
13.6 X 10-3 
6.0 X 10-3 


«.U_o42 

8 7+0.42 
''■'-0.42 
8 S+O ™ 
°-^-0.48 

8 fi+o-se 

S.t)_o35 

7 20+'''36 

'■^"-0.36 



Table 2. X-ray parameters for the five clusters used in the combined Ho measurement. 
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Figure 1. The full NORAS X-ray cluster sample and the region (white) of 
our SZ subsample. The lines coiTespond to flux limits of 1, 3 and 5 x 10^'~ 
erg s^' cm^^, which are the approximate lower detection limit, the 90% 
completeness limit, and the limit of the SZ sample respectively. 
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Figure 3. Source-subtracted visibilities (real part only) for (left) A697, (centre) A1914 and (right) A2218 (filled points with error bars), binned azimuthally, 
compared with the best-fitting model in each case (open points joined by lines). The lack of smoothness in the models is due to the inclusion of points from 
different baselines with different model values (due to eUipticity) in the same bin. 
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Figure 4. (Left)CLEANed map of all the radio data on A697 (contours) superposed on the ROSAT HRI image (greyscale). The contour interval is 
60/jJybeam~^; dashed contours are negative. The CLEAN restoring beam of 101" x 61" is shown (inset). (Right)CLEANed map of the source-subtracted 
radio data on A697, with a gaussian taper of width 2000X applied (contours) superposed on the ROSAT HRI image (greyscale). The contour interval is 
60;uJybeam^' ; dashed contours are negative. The CLEAN restoring beam of 76" x 39" is shown (inset). 
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Figure 5. (left) Full-resolution radio image of A1914 (contours) overlaid on a VLA 1.4 GHz image (greyscale). The contours are at —3,-2,2,3,4,5,6,7,8 x 
35/Jybeam~',; the beamsize of 69 x 33 arcsec, PA. 5° is shown (inset). The steep-spectrum halo source can be seen in the VLA image just to the east of 
the cluster centre - no emission is apparent in the RT image however, where the SZ decrement can also be seen, (right) Source-subtracted CLEANed low- 
resolution radio image of A1914 showing the SZ decrement, overlaid on the ROSAT PSPC image (greyscale). The contour interval is SO^uJybeam^', dashed 
contours are negative. The beam (inset) is 175" x 108", PA 4°. 
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Figure 6. (left) Full-resolution radio image of A2218 (contours) overlaid on the NVSS 1.4 GHz image (greyscale). The contours are at 
—4,— 3,— 2,2,4,8,16,32 x 35/jJybeam~'; the beamsize of 41 x 36 arcsec, PA —9° is shown (inset). The SZ decrement can be clearly seen even in this 
high-resolution image, (right) Source-subtracted CLEANed low-resolution radio image of A2218 showing the SZ decrement, overlaid on the ROSAT PSPC 
image (greyscale). The contours interval is 80// Jybeam"^ , dashed contours are negative. The beam (inset) is 130" x 120", PA —5°. 
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